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Abstract 


Transient  acoustic  signals  are  studied  in  compound  waveguides 
consisting  of  a  resonant  cavity  attached  to  a  semi-infinite  cylindrical 
waveguide.  The  signals  are  shown  to  have  the  asymptotic  form 

00 

v(t,x,y)  -  l  v.  (t,y)$.  (x)  ,  t  -*■  “»  , 
j-0  J  J 

where  x  ■  (x  ,x  )  are  coordinates  in  the  cylinder  cross-section,  y 
1  2 

is  a  coordinate  along  the  cylinder  and  t  is  a  time  coordinate.  Here 
<|>j  (x)  is  an  eigenfunction  for  the  cylinder  cross-section,  with  eigen¬ 
value  y^  ,  and 

Vj(t,y)  -  9(t,y,yj)  F^  (y^  y/(t2  -  y2)  ) 

where  9(t,y,y)  is  a  universal  factor  and  Cp)  characterizes  the 
momentum  distribution  of  mode  j  .  It  is  shown  that  if  both  the  signal 
sources  and  observation  point  are  far  from  the  resonator  then 

Fj (P)  -  Fjir  (p)  +  FjC  (p) 

is  the  direct  wave  that  would  exist  if  no  resonator  were 

-  SFj1'  <p)  -  l  SJm(p)F"r  </pJ+V“„)  • 

p2+y  -y  >0 
m  n 

A 

S  is  the  S-matrix  for  the  compound  waveguide  and  may  be  calculated 
from  the  model  functions. 


where  F 


,dir 

j 


present  and 


FjC  (P) 


Figure  1.  Cylindrical  waveguide  coupled  to  a  resonator. 

The  walls  of  the  waveguide  are  assumed  to  be  rigid.  The  sources 
of  the  transient  sound  fields,  or  "signals,"  are  assumed  to  be  localized 
in  a  bounded  portion  of  the  waveguide  and  to  act  for  a  finite  interval 
of  time.  The  goal  of  the  work  is  to  calculate  such  acoustic  signals  and 
to  analyze  how  their  structure  depends  on  the  sources  and  the  geometry 
of  the  waveguide.  Particular  attention  is  given  to  the  cases  in  which 
the  signal  sources  or  observation  point,  or  both,  lie  in  the  cylinder 
and  are  far  from  the  resonator. 

It  is  well  known  that  in  the  cylindrical  portion  of  the  waveguide 

each  signal  can  be  decomposed  into  a  series  of  modal  waves.  These  waves 

o 

are  calculated  below  and  are  shown  to  be  asymptotically  Independent  for 
large  times.  Moreover,  the  form  of  the  modal  waves  is  shown  to  be 


determined  by  the  geometry  of  the  cylinder.  Only  the  fine  structure  of 
the  modal  wave  profiles  varies  with  the  sources  and  the  geometry  of  the 
resonator. 

The  theory  of  the  waveguides  of  Figure  1  will  be  developed  by 
perturbation  theory,  beginning  with  the  special  case  of  the  simple 
waveguide  consisting  of  the  semi-infinite  cylinder,  without  resonator, 
terminated  by  a  plane  cap.  The  general  case  of  a  cylinder  plus 
resonator,  depicted  in  Figurel,  will  be  called  a  compound  waveguide. 

Acoustic  signal  propagation  in  both  simple  and  compound 
waveguides  will  be  analyzed  by  means  of  normal  mode  expansions.  For 
simple  waveguides  the  normal  modes  can  be  constructed  explicitly  by 
separation  of  variables.  For  compound  waveguides  they  are  constructed 
by  a  perturbation  method  based  on  those  for  simple  waveguides.  The 
asymptotic  wave  functions  that  describe  the  signals  for  large  times  are 
calculated  from  the  normal  mode  expansions. 

The  asymptotic  wave  functions  for  a  waveguide  are  characterized 
by  a  sequence  of  functions  that,  physically,  describe  the  momentum 
distributions  of  the  modal  waves.  With  each  waveguide  is  associated  a 
scattering  operator,  or  S-matrix,  that  operates  in  the  space  of  these 
momentum  distributions.  The  final  result  of  this  paper  reveals  the 
significance  of  the  S-matrix  in  the  analysis  of  acoustic  signal  struc¬ 
ture  in  waveguides.  It  states  that  the  momentum  distribution  of  the 
signal  generated  in  a  compound  waveguide  by  sources  far  from  the 
resonator  is  simply  the  image  under  the  S-matrix  of  the  momentum 
distribution  of  the  signal  generated  by  the  same  sources  in  the 
corresponding  simple  waveguide. 


Normal  mode  expansions  for  waveguides  and  the  associated  S-matrlx 
were  developed  by  C.  Goldstein  In  the  period  1969-74;  see  [1,2, 3,4]*  In 
1975  Goldstein's  results  were  extended  by  W.  C.  Lyford  [5],  using  results 
of  the  abstract  theory  of  scattering.  Lyford  also  presented  results  on 
asymptotic  wave  functions  for  waveguides  in  [6].  In  1977  the  author 
presented  in  [11]  an  exposition  of  the  theory  of  normal  mode  expansions 
and  asymptotic  wave  functions  for  the  more  general  case  of  several 
semi- inf inite  cylindrical  waveguides  coupled  by  a  resonator;  see 
Figure  2. 

The  purpose  of  this  paper  is  to  present  a  new  construction  of 
the  S-matrix,  based  on  the  theory  of  asymptotic  wave  functions,  and  to 
apply  the  S-matrlx  to  the  construction  of  asymptotic  wave  functions  due 


Numbers  in  square  brackets  indicate  references  from  the  list  at  the 
end  of  the  paper. 


to  sources  remote  from  the  resonator.  The  work  is  based  on  the  earlier 
literature,  primarily  the  exposition  [11],  and  the  author's  monograph  on 
the  closely  related  theory  of  scattering  by  diffraction  gratings  [12] . 

It  will  be  seen  that  the  results  presented  below  can  also  be  derived  for 
acoustically  soft  boundaries  (Dirichlet  condition)  and  elastic  boundaries 
(Robin  condition) ,  and  for  the  general  case  of  a  resonator  with  several 
cylinders.  However,  to  simplify  the  exposition  only  the  case  of  a  single 
cylinder  and  rigid  boundary  is  treated  here. 

The  remainder  of  the  paper  is  organized  as  follows.  §2  presents 
a  formulation  of  the  propagation  problem  for  acoustic  signals  in 
compound  waveguides  as  an  initial -boundary  value  problem  for  the  wave 
equation,  together  with  its  solution  by  a  simple  Hilbert  space  method. 

S3  develops  the  normal  mode  expansion  for  simple  waveguides.  In  §4 
the  normal  mode  functions  for  compound  waveguides  are  defined  and  the 
normal  mode  expansions  for  this  case  are  presented.  In  §5  the  normal 
mode  expansions  of  acoustic  signals  in  waveguides  are  derived  and  their 
asymptotic  wave  functions,  for  large  times,  are  calculated.  Sections  2 
through  5  present  a  review  of  concepts  and  results  from  [11].  The 
exposition  is  therefore  concise  and  without  proofs.  The  new  results  of 
this  paper  are  contained  in  §6  and  §7.  In  §6  the  S-matrix  of  a  compound 
waveguide  is  defined  and  then  constructed  by  means  of  the  asymptotic 
wave  functions  of  §5.  In  §7  the  acoustic  signals  generated  by 
prescribed  sources  in  simple  and  compound  waveguides  are  compared.  The 
principal  result  of  this  section  states  that,  for  sources  far  from  the 
resonator,  the  momentum  distribution  of  the  scattered  signal  is  simply 
the  image  under  the  S-matrix  of  the  momentum  distribution  generated  by 
the  same  sources  In  the  simple  waveguide. 


.  The  Propagation  Problem  and  Its  Solution. 


The  notation  that  is  used  in  the  remainder  of  the  paper  is  fixed 
inthis  section.  Acoustic  fields  in  waveguides  will  be  described  by  real 
valued  acoustic  pressure  functions  u(t,X)  where  t  €  R  is  a  time  coordinate 
and  X  -  (x,y)  -  (x^Xj.y)  e  R3  represents  a  triple  of  rectangular 
coordinates  in  space.  The  coordinate  axes  are  assumed  to  be  chosen  in 
such  a  way  that  the  simple  waveguide  occupies  the  semi-infinite  cylinder 

(2.1)  ft0  *  G  x  R0  -  {X  :  x  €  G  and  y  >  0} 

where  G  is  a  bounded  domain  in  the  ,x2 -plane  and  Rj  •  (y  :  y  >  0}. 

The  corresponding  compound  waveguide  domains  are 

(2.2)  ft  -  ft„  U  K 

where  K  is  a  bounded  domain  and  ft  is  connected,  and  hence  is  a  domain. 

The  boundary  of  ft  is  denoted  by  9ft.  It  must  be  mentioned  that  the  local 
structure  of  9ft  is  not  completely  arbitrary  since  ft  will  be  required  to 
have  the  local  compactness  property  of  [11,  p.  408].  A  simple  geometrical 
property  that  is  sufficient  to  guarantee  this  is  the  finite  tiling 
property  of  [10] .  All  the  simple  piecewise  smooth  boundaries  that  arise 
in  applications,  such  as  unions  of  polyhedra,  and  finite  sections  of 
cylinders,  cones,  spheres  and  disks,  may  be  shown  to  have  this  property. 

The  acoustic  pressure  u(t,X)  is  the  solution  of  an  initial¬ 
boundary  problem  for  the  wave  equation  which,  in  its  classical  formu¬ 
lation  ,  reads 


0  for  t  >  0,  X  e  ft, 


(2.3) 


(2.5)  u(0,X)  «  f(X)  and  — U-^?^  ■  g(X)  for  X  €  ft 

where  V  and  A  are  the  gradient  and  Laplace  operators  in  R3 ,  respectively, 
while  v  is  a  unit  normal  field  on  3ft.  The  functions  f(X)  and  g(X)  in 

(2.5)  characterize  the  initial  state  of  the  acoustic  field.  They  are 
assumed  to  be  given  or  calculated  from  the  prescribed  wave  sources  in  ft. 

A  general  theory  of  the  initial-boundary  value  problem  (2.3)  - 

(2.5) ,  guaranteeing  the  existence  and  uniqueness  of  the  solution  for 
arbitrary  domains ,  was  given  in  [9].  Of  course,  for  arbitrary  domains 
the  boundary  condition  (2.4)  must  be  understood  in  a  generalized  sense. 
Here  a  simple  alternative  to  the  method  of  [9]  will  be  based  on  the 
acoustic  propagator  A  in  the  Hilbert  space  J€  *  L2 (ft) .  A  is  the 
selfadjoint  realization  of  the  negative  Laplacian  defined  by 

(2.6)  D(A)  »  3f  n  {u  :  Vu  and  Au  are  in  3C,  *  0  on  3ft} 

and 

(2.7)  Au  *  -Au  for  all  u  €  D(A) . 

The  differential  operators  V  and  A  in  (2.6)  are  to  be  interpreted  in  the 
sense  of  distribution  theory.  Of  course,  if  3ft  is  not  smooth  then  the 
boundary  condition  in  (2.6)  is  to  be  interpreted  in  the  generalized  sense 
of  [9].  With  these  conventions  it  can  be  shown  that 

(2.8)  A*  =  A  >  0  in  JC; 


see  [10,11]  for  details. 
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With  the  above  definitions  the  problem  (2.3)  -  (2.5)  may  be 
reformulated  as  an  initial  value  problem  for  a  function  t  -*■  u(t,*)  €  3C; 
namely 

(2.9)  +  Au  *  0  for  t  >  0 

(2.10)  u(0)  -  f  and  .  g. 

A  formal  solution  is  given  by 

(2.11)  u(t,*)  ■  (cos  t  A1^2)  f  +  (A-1^2  sin  t  Al/,z)  g. 

The  spectral  theorem  Implies  that  (2.11)  defines  the  "solution  in  X"  of 
(2.9),  (2.10)  for  all  f,g  €  X.  It  is  also  the  unique  "solution  with 
finite  energy"  of  [9]  for  all  initial  states  (f,g)  with  finite  energy. 
This  is  equivalent  to  f  €  D(A1^2),  g  €  X. 

It  will  be  convenient  to  represent  the  solution  u(t,X)  as 

(2.12)  u(t,X)  -  Re  {v(t ,X) } 

where  v(t,X)  is  the  complex-valued  function  defined  by 

-itA1/2 

(2.13)  v(t , •)  -  e  itA  h 
and 

(2.14)  h  -  f  +  i  A-I/z  g 

This  representation  is  valid  if  f  and  g  are  real  valued  and  f ,  g  and 
A  ^2  g  are  in  X.  It  is  sufficient  to  consider  such  initial  states  since 
the  spectral  theorem  implies  that  they  are  dense  in  X. 

4’ 


.  Normal  Mode  Expansions  for  Simple  Waveguides. 

The  acoustic  propagator  for  the  simple  waveguide  £)<>  *  G  x  &o 
will  be  denoted  by  A<j .  It  is  selfadjoint  and  non-negative  in  the  Hilbert 
space  X0  ■  L2(fi0).  The  spectral  family  of  Ag  was  calculated  in  [11]  and 
shown  to  be  absolutely  continuous.  Here  the  generalized  eigenfunction, 
or  normal  mode,  expansion  for  A0  is  reviewed  briefly. 

The  normal  mode  functions  for  Aq  may  be  found  by  separation  of 
variables  [11,  p.  418ff] .  They  are  defined  by 


(3.1)  i/jj (x,y ,p)  ■  I— I  <f)j(x)  cos  py,  p  >  0,  j  -  0,  1,  2,***  , 


where  the  functions 


(3.2) 


d>0  (x)  -  |G|~1/2  ,  <h(x),  (x)  ,  * #  • 


are  the  eigenfunctions  of  two-dimensional  negative  Laplacian 


AG  ■  -3z/3x2  -  32/3x2  in  L2(G),  with  Neumann  boundary  condition.  The 

N 

corresponding  eigenvalues  of  A„,  repeated  according  to  their 

V* 

multiplicity  will  be  denoted  by 


(3.3) 


Uo  -  0  c  Uj  <  U2  <  •*' 


It  is  known  that  each  y^  has  finite  multiplicity  and  y^  -*•  00  when  j  -*■  <*>. 
Moreover,  the  smallest  eigenvalue  is  simple  and  <j>0  (x)  *  const.  | G j  is 
the  2-dimensional  Lebesgue  measure  of  G.  The  eigenvalue  problem  for  the 
in  classical  form,  is 


(3.4) 


(A  +  wj(p))  t|>°(X,p)  -  0,  X  6  fl0. 
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(3.5)  3t{)®(X,p)/3v  -  0,  X  6  3a0, 
where 

(3.6)  tUj  (p)  -  (p2  +  yj)l/2  >  uj/2  J-0,1,  2,-  . 

The  normal  mode  expansion  theorem  for  Ao  states  that  every 
h0  S  Kq  can  be  written  as 

00 

(3.7)  h0(x,y)  -  l.i.m.  £  h  .(y)  <M*) 

j-0  2  2 

where  l.i.m.  denote  convergence  in  3C0  [11].  The  functions 
h0j  e  1*2  (Ro )  are  defined  by 

hV12 

(3.8)  hQ  j  (y)  =  l.i.m.  jo  cos  py  h„j(p)  dP  » 

where 

(3.9)  h.  ,(P)  -  l.i.m.  I  i|;°(X,p)  h„(X)  dX  . 

3  JB0  3 

In  (3.8)  and  (3.9),  l.i.m.  denotes  convergence  in  L2  (I>0).  The  Parseval 
relation  for  the  expansion  (3.7)  -  (3.9)  is  contained  in  the  theorem 
that  the  linear  mapping  $  defined  by 

(3.10)  4>h0  *  {h00  ,h01  ,h02,***} 
defines  a  unitary  operator 


The  last  notation  means  that  K0  is  the  direct  sum  of  a  countable 
sequence  of  copies  of  L2(R0). 


The  representation  $  is  a  spectral  mapping  for  A0 .  More 
precisely,  one  has 


(3.12) 


v0 (t,X) 


e-itA^2 


h0  (X) 


l.i.m.  I  v  ,  (t,y)  <|>,(x)  in  JC( 
j-0  02  2 


where 


(3.13)  vflj(t,y) 


l.i.m. 


fCO 

o 


cos  p  y  e 


-itoij  (p) 


fioj(p)  dp 


Normal  Mode  Expansions  for  Compound  Waveguides. 


For  compound  waveguides  Q  •  K  U  the  propagator  A  has  two 
distinct  families  of  normal  mode  functions.  They  are  denoted  here  by 
{4>j(X,p)  :  p  >  0  and  j  *0,  1,  2,»»»)  and  {^(X.p)  :  p  >  0  and 
j  ■  0,  1,  2,***}.  Both  families  satisfy  the  differential  equation  and 
boundary  condition: 


(4.1)  (A  +  wj(p))  ijjJ(X,p)  ■  0  for  X  €  Q  , 

(4.2)  3\Jk (X,p) /3v  -  0  for  X  €  . 


The  families  are  distinguished  by  the  condition  that  in  (20  the  fields 


,±sc 


(4.3)  (X,p)  -  l|>J(X,p)  -  lp](X,p),  X  6 


represent  outgoing  waves  for  and  incoming  waves  for  .  This  is 


+sc 

defined  by  the  condition  that  the  developments  of  \)>,  in  the  transverse 


j 


eigenfunctions  (<j>m(x)}  have  the  forms 


^sc(x,p) 


(4.4) 


1  -  ±  ±iy^(p)-ym 

JW>u  Vrt  * 

j  m 


i  P  +  Vv^(p) 

+  -jLsv  T5"<p)  e  ' 


The  two  summations  in  (4.4)  are  over  the  sets  of  integers  m  >  0  for  which 
the  indicated  inequalities  hold.  Note  that  the  first  sum  is  finite 
because  u  -*■  »  when  m  ■+•  «>.  The  coefficients  {TT  (p)  :  co?(p)  >  u  }  will 

m  jm  j  m 

be  called  the  scattering  amplitudes.  It  will  be  shown  that  they 

determine  the  S-matrix  for  the  compound  waveguide. 
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F.  Rellich  [8]  was  the  first  to  show  that  for  compound  waveguides 
ft  *  K  'J  ft0  the  acoustic  propagator  A  may  have  point  spectrum 
(\n)  :  1  -  n  <  **  5  +00}  and  eigenfunctions 

(4.5)  ¥(n)  e  D(A)  C  X 

such  that 


(4.6) 


A^(n)  *  X(n)  *(n) 


and 


(4. 7)  ^(n)(x>  * 

The  wave  functions 
acoustic  fields 


U  >X.  . 
m  (n) 


-y/u  -X.  . 

1  m  (n)  .  ,  . 
c  e  <p  (x)  . 

nm  rm 


may  be  called  trapping  modes. 


The  corresponding 


M  -itX1/2. 

(4.8)  v(t,X)  -  l  c^  e  n)  ip(n)  (X) 


represent  standing  waves  in  ,  by  (4.7).  In  particular,  the  energy  of 
the  wave  function  (4.8)  does  not  propagate  to  »  and  hence  plays  no  role 
in  the  scattering  theory  for  A.  In  the  remainder  of  this  paper,  to 
simplify  the  notation,  it  is  assumed  that  A  has  no  trapping  modes.  In 
the  general  case  the  results  derived  below  hold  in  the  space  of  states 
orthogonal  to  the  trapping  modes. 

For  propagators  A  with  no  trapping  modes  the  normal  mode 
expansion  theorem  states  that  every  h  6  Jf  can  be  written  as 


(4.9)  h(X)  =  l.i.m.  \  h*(X) 

j-0  J 


>  -  a  -  i  A 


15 


where  the  series  converges  in  X.  The  components  h“  £  3C  are  defined  by 


(4.10) 


h“(X)  -  1.1 


.m.  i^(X,p)  h*(p)  dp  , 


with  convergence  in  X,  where 


(4.11) 


hT(p)  -  1.1 


m.  f  C(X, 


p)  h(X)  dX 


in  L2(R0).  The  orthogonality  and  completeness  of  the  expansions  are 
expressed  by  the  theorem  that  the  linear  operators  $+  and  $>  defined  by 


(4.12) 


are  unitary  operators  from  X  to  JC. .  The  spectral  property  of  and  $ 


is  characterized  by  the  representation 


(4.13)  v(t,X)  -  e'itA‘'‘ h(X)  *  l.i.m.  J  v*(t,X)  . 


1/2 


j-0 


where 


(4.14)  v“(t,X)  -  l.i.m. 


*  +  -itd).  (p)  A+ 

^ . (X,p)  e  J  h  (p)  dp 

0  J  J 


The  convergence  in  both  (4,13)  and  (4,14)  is  in  Jf. 


.  Asymptotic  Wave  Functions  for  Waveguides 


The  theory  of  asymptotic  wave  functions  for  waveguides,  as 
developed  in  [6]  and  [11],  is  reviewed  briefly  in  this  section. 

Consider  first  the  wave  function  v0(t,X)  for  the  simple  wave¬ 
guide  defined  by  (3.12),  (3.13).  The  modal  wave  voj  can  be  written 


(5.1) 


voj(t,y)  -  v°Jfc(t,y)  +  vjj(t.y) 


where 


(5.2)  v^(t,y)  -  v£(t,-y)  -  ^ 


i(yp-tu)j  (p)) 


fio;)(p)dp  • 


Here  and  in  what  follows  the  l.i.m.  notation  of  (3.12),  (3.13)  is 


omitted,  for  brevity,  but  is  always  to  be  understood. 

The  wave  functions  of  v0^1"  and  v*?  have  the  form 

o  j  oj 


(5.3)  w  (t,y,y,h)  « 


ei(±yp-cw(p,u))  h(p)  Jp 


where  h  €  L2(R„)  and  uj(p,y)  “  (p2  +  y)1^2  •  For  y  >  0,  (5.3)  defines  a 
dispersive  wave  with  group  speed 


(5.4) 


U(p,y) 


3to(p,y)  _  p 

3p  ~  (pz4y) 


In  [11]  the  method  of  stationary  phase  was  applied  to  w+.  The  resulting 


estimate  is 


w  (t,y,y,h) 


(5.5) 


i(yp-tu>(p,y)-ir/4) 

n  « _ 

tJ  (t3u(p,u)/3p)1/2 


0(t,y,y)  h j 


U  -yzJ 


2  )  1/2 


where 


(5.6) 


Characteristic  function  of  {(t,y)  :  0  <  -2-  <  1} 


(5.7)  0(t,y,y)  ■  X  tu[u(t2  -  y2)  J-^1*  e-i[w(t2-y2)  1/2  +ir/4] 


For  u  *  0,  (5.3)  reduces  to 


(5.8)  w±(t,y,0,h)  -  (2irj-i/a  r  ei(±y_t)p  h(p)  dp  -  F'1  h(±y-t) 


where  F  is  the  usual  Fourier  transform.  In  this  case  the  waves 
non-dispersive  and  it  is  appropriate  to  define 


(5.9) 


w  (t,y,0,h)  -  F~l  h(y  -  t)  . 


With  these  definitions  the  results  of  [11]  imply  that,  for  all  h  e  L2(RQ) 
and  all  u  >  0  one  has 


(5.10)  lim  I w+(t, •  ,u,h)  -  w°°(t, •  ,y ,h)l  ,  .  -  0  , 


(5.11) 


lim  I w”(t, • ,u,h)!  ,  -  0 

t-Ko  L2'K0' 


These  results  are  directly  applicable  to  the  partial  waves  v0^ 

defined  by  (5.1),  (5.2).  It  follows  that  v*"(t,*)  -►  0  and 
Out  oo 

(t,*)  -  v0j(t,*)  -*■  0  in  L2(Ro)  when  t  00  where 

y  y 2  y 

(5.3.2)  v® ,(t,y)  -  0(t,y,y.)  F  —  j  ---  ,  j  -  1,  2,***, 

2  2  [(t2-y2) 1/2 

(5.13)  v®0 (t»y)  ■  (f-1  F80)(y-t)  , 
and 

(5.14)  Foj(p)  -  h0j(p)  for  j  *  0,  1,  2 , • • *  . 

Moreover,  it  was  shown  in  [11]  that  if 

OO 

(5.15)  v<o (t,x)  -  l  v“  (t,y)  4>.(x) 

j-0  J  J 

then  the  series  converges  in  3C0  and 

(5.16)  lim  Iv„(t,*)  -  v“(t,*)HjC  *  0  . 

t**» 

OO 

It  is  clear  from  these  results  that  the  asymptotic  wave  function  vQ  is 
characterized  by  the  sequence  {F„j}  of  modal  wave  momentum  distributions 
or,  equivalently,  the  element 

(5-17)  F0  -  {F0j}  =  {h0.}  =  $h0  €  jf0  . 

The  same  analysis  may  be  applied  to  the  wave  function  v(t,X)  for 
the  compound  waveguide  defined  by  (4.13),  (4.14).  If  the  "incoming" 
normal  modes  {$.}  are  chosen  the  partial  waves  have  the  form 


(5.18)  v“(t,X)  -  jB  ^(X.p)  e  J  h^p)  dp  . 


For  X  €  J)0  the  decomposition  (4.3),  (4.4)  for  may  be  written 


VX*P)  -  (2^p7r-eiPy  Vx) 


(5.19) 


1  “  -iy/a)f(p)-u 

*  Jn  SJ“(P)  6  *m(x)  VP) 


1  “  -y/u  -wif(p) 

+  (2»j-i7*-  JQ  Sjm(P>  e  Vx)  Xj'm(p) 


where 


(5.20) 


s>> '  5j»  +  T>> 


(5.21)  X.,  (p)  ■  Characteristic  function  of  (p  :  oo*(p)  >  y  } 


and  Xjm(p)  +  Xjm(p)  -  1.  Substituting  (5.19)  into  (5.18)  gives 


(5.22) 


.  ,;.<■«  +  v?-  + 


where  the  three  terms  on  the  right  come  from  the  corresponding  terms  of 


(5.19) 


Clearly 


(5  *23)  Vj’out(t,X)  -  w+(t,y,y..  ,h~)  ^(x)  . 


Moreover. 


m»0 


1 

(2ir) 


n  S.(p) 
o  j® 


-iy/o)^ (p)-u ~  itu. (p)  ^ 
e  3  m  3  hj(p) 


Xjm<P)<>P 


rm 


On  making  Che  change  of  variable 

(5.25)  P*2  -  OJj(p)  -  Um  *  P2  +  Wj  -  Hm 
In  (5.24)  one  finds  ChaC 

00 

(5.26)  v“’in(t,X)  -  l  w‘(t,y,Mm,h"  )  *  (x) 

J  m-0  J 

where  the  function  h.  is  defined  by 

jrn 

(5.27)  VP'>  "  Xjn<P)  STm(p)  h'(p) 

and  p  and  p'  are  related  by  (5.25).  It  is  clear  from  (5.11)  and  (5.26) 
that,  formally, 

(5.28)  lim  v“’in(t,*)  -  0  in  Jf0. 

Rigorous  proofs  of  this,  and  the  related  property 

(5.29)  lim  v"’a(t,0  =  0  in  JC0  , 

were  given  for  the  closely  related  case  of  diffraction  gratings  in 
[12,  Part  2].  (5.23),  (5.28)  and  (5.29)  imply  that 

vT(t,*)  -  vT,C°(t,-)  -  0  in  X 


(x) 


(5.30) 


_ _  I  II  I  .  .  II  J  I  .  I  I  I  1  -1 -  —  - - -y  —  ---  vy  «Y,1 

.*.  ■' .  »*.  ■  .  -  .  *'V*  «_  y _■  .-.  ,  -  _ - ,  -.-.•»  •_.  v» >y.  ■  i  ■  .  vJ 
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(5.31) 


v”,0°(t,X)  -  x(X)  w“(t,y,Uj,h”)  <J>j(x) 


and  X  Is  the  characteristic  function  of  H0 .  Moreover,  if 


(5.32) 


v-’“(t,X)  -  I  vT’“(t,X) 
j-0  3 


then  it  follows  as  in  [11]  that  the  series  converges  in  3 C  and 


(5.33) 


lim  lv(t,*)  -  v  *  (t,*)ljf  *  0 


t-»oo 


for  every  h  €  X.  Note  that 


(5.34)  Vj’  (t,X)  -6(t,y,UJ)  Fj 


l(t2-y!)1/s  J 


4>j(x)  ,  j  -  1,  2,  3,***, 


(5.35) 


v7”(t,X)  -  (F-1  Fs)(y-t)  |g|~1/2 


where 


(5.36) 


Fj(p)  -  h'(p)  for  j  *  0,  1,  2 , • • •  . 


QQ 

Ulus  v  *  is  characterized  by  the  sequence  { Fj }  of  modal  momentum 


distributions  where 


(5.37) 


F  -  {Fj}  =  {h~}  *  <Th  e  ic0 


The  convergence  of  v(t,*)  to  zero  in  L2 (K) ,  which  is  implied  by  (5.33), 
follows  from  the  local  compactness  property  of  [11,  p.  408];  see 
[11,  Theorem  6.16]. 


A 
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The  results  (5.16)  and  (5.33)  imply  that  each  wave  function 
v(t,*)  for  the  compound  waveguide  is  asymptotically  equal  in  for 
t  -*■  00 ,  to  a  sum  of  modal  waves  for  the  simple  waveguide. 


k 


^  >t  . 


.  ii  mi i  ■m. ■  r  .■ 


6.  The  S-Matrix. 


(6.1) 


The  wave  operators  W+  ■  W+(A^2  .A1^2  ,J)  :  JC  ■+■  3f#  are  defined  by 

„  ..  itAo1/2  T  -itA1/2 

W+  -  s-lim  e  J  e 

—  f^i00 


where  J  :  X  ■*  X0  is  the  bounded  operator  defined  by  (Jh)(X)  ■  h(X)  L 
and  s-lim  denotes  strong  convergence.  The  existence  of  W+  follows  from 
(5.16),  (5.17),  (5.33)  and  (5.37).  More  precisely,  these  results  imply 
that 


(6.2) 


lim  D  e 
t-*» 


itA°1/2  J  e~itAl/‘  h  -  h.l 


1/2 


0  5f„ 


provided  F„  »  4>h„  *  $  h  «  F.  This  implies  the  representation 


(6.3)  W, 


The  analogous  results  for  t 


(6.4) 


give 


W  -  $*  $+  . 


These  results  and  the  unitarity  of  $  and  4>~  imply  that  W+  are  unitary 
operators  from  X  to  XQ . 

The  scattering  operator  S  for  the  pair  is  the  unitary  operator 
in  X0  defined  by 


(6.5) 


s  -  w+  w_  . 


Combining  (6.3),  (6.4)  and  (6.5)  gives 


(6.6) 


*  ^ 

S  -  <t>  S  4> 
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A  /\ 

where  S  is  the  unitary  operator  in  the  direct  sum  space  3C0  of  (3.11) 
that  is  defined  by 

(6.7) 

A 

S  will  be  called  the  S-matrix  for  A  (or  the  the  waveguide  Q) . 

The  purpose  of  this  section  is  to  calculate  an  explicit 

representation  of  S  based  on  the  scattering  amolitudes  {S,  (p)}  of  the 

Jm 

normal  mode  functions  ^(X,p),  defined  by  (4.4)  and  (5.20).  To  this 
end  let  h  €  Jf  and  define  elements  h+  and  h  in  3C0  by  (4.12);  i.e., 

a+  a+  a+  _  + 

(6.8)  h-  =  {hj ,hj ,hj ,*••}  =  h  . 

Then  it  follows  from  the  unitarity  of  $+  and  that 

(6.9)  h  *  ($  $+*)  $+  h  *  S  h+  . 

This  relation  will  be  used  to  calculate  S.  To  see  how  this  can  be  done 

recall  that  the  calculation  of  the  asymptotic  wave  function  for 
-itA1/2 

v(t,*)  *  e  h  was  based  on  the  incoming  representation  (5.18)  and 

gave  the  representation  (5.37),  i.e.,  F  =  h  ,  for  the  asymptotic 
momentum  distribution.  A  second  calculation  of  F,  based  on  the  outgoing 
representation,  will  now  be  shown  to  yield  a  representation  of  §. 

The  outgoing  representation  of  v(t,*)  is 


(6.10) 


where 


v(t,X)  -  l  v+(t,X) 

J-0  J 


v+(t,X) 


fCO 

<(X,P) 
0  J 


-itu),  (p) 

;  J  b^(p)  dp 


(6.11) 


For  X  S  U0  the  decomposition  (4.3),  (4.4)  for  ^  may  be  written 


(6.12) 


where 


♦Jtt.p)  -  J2W177'  e"iPy  *j(x) 


1  r  +  iy^(p)-u_ 

+  ('^f72'  J0  Sjm(p)  6  *m(x)  Xjm(P) 

1  n  +  -y/y-u^  (p) 

+  72^  i  Sjm(p)  e  Vx)  Xjm(P> 


(6.13) 


S]m(p)  ’  5jm  +  T]>> 


Substituting  (6.12)  into  (6.11)  gives 


(6.14) 


+  +,in  +,out  .  +,a 

Vj-Vj  +Vj-  +Vj  , 


where  the  three  terms  on  the  right  come  from  the  corresponding  terms  of 
(6.12).  The  term 


(6.15) 


vj’in(t,X)  *  w-(t,y,y;j,h|)  $j(x) 


tends  to  zero  in  3C0  when  t  -*■  00 ,  by  (5.11).  Moreover,  the  method  of 
[12,  Part  2,  §14]  may  be  used  to  show  that  v*,a(t,*)  0  in  3f0  when  t  » 

To  calculate  the  behavior  for  t  -►  00  of  v(t,X)  from  (6.10)  -  (6.14 
it  is  convenient  to  begin  with  the  special  case  of  an  h  e  Jf  such  that 


(6.16) 


hj(p)  =  g(p) ,  n  fixed,  and 


(6.17) 


supp  g  C  [0,M],  M  fixed  . 


For  such  h 


(6.18)  v(t,X)  -  v+(t,X)  -  v+’out(t,X)  +  0(1) 

where  0(1)  -*-0  in  3C0  when  t  -*■  «.  Moreover, 


(6.19) 


supp  Xqjjj  ■  ['Max.  (0,um-yn)  ,®) 


and  hence 


v+'°ut(t,X) 

(6.20) 


U<y+M2 
m  n 


TJ7- 


M  iy/w2 (p)-y  -  itu(p) 

Sm/P)  e  (j>m(x)g(p)dp 

/Max(0,y  -y) 

m  n 


The  change  of  variable  p  -*■  p'  where 


(6.21) 


wn(p)  -  p2  +  un  •  p'2  +  ym  -  m2(p’) 


gives  the  representation 


v+’°Ut(t,X) 


(6.22) 


£ 


(2rr)  ^ 


-y 

+  i(yp'-ta,m(p*)) 

Snm(p)  6  <l>m(x)  g(p)  dp' 


y  <yn4M" 

m—  n 


/Max(0,y-y) 
n  m 


where  in  the  last  integrand  p  »  /p ' 2  +  yffl  -  yn .  The  definition  (5.3)  of 
w+  implies  that  (6.22)  can  be  written 


(6.23) 


v+’°Ut(t,X)  = 


w+(t,y,y .h*)  4>  (x) 
2  m  nm  m 


where 

(6.24)  .  S*  (p)  g(p>  (P'>  • 


Note  that  the  lower  limit  of  Integration  In  (6.22),  together  with  (6.19), 
imply  the  correctness  of  the  factor  ^(p')  In  (6.24). 

Equations  (6.18),  (6.23)  and  (5.10)  imply  that  if 


(6.25) 


v+ “(t,X)  -  I  „  wa6Ct,y,um,h^m)  <f>m(x) 


U<U  4M2 
m-  n 


then 


(6.26) 


lim  lv(t,*)  -  v+,°°(t,*)i 

t-wo 


Combining  this  with  (5.33)  gives 


(6.27) 

where 


Iv‘*“(t,-)  -  v+’“(t,.) 


'K 


(6.28) 


v-,0°(t,X)  -  l  w+(t,y,um,h~)  d>m(ac) 
m=0 


Equations  (6.25),  (6.27)  and  (6.28)  imply  that 


(6.29) 


h  =  h+  for  m  *  0,  1,  2,  ■ 
m  nm 


This  is  a  consequence  of  a  simple  lemma  which  states  that  if 

h  =  {h  }  €  and  if 
m  u 


(6.30) 


w(t,X)  -  \  w°°(t,y,ytn,hra)  4>n(x)  -*•  0 

m-0 


in  when  t  -*■  °»  then  h  *  0  in  JC0  ;  i.e.,  0  ■  h0  *  hj  =  h2  «  •••  .  To 


verify  the  lemma  note  that  (6.30)  is  equivalent  to  the  conditions 


(6.31) 


w  (t,y,ym,hm)  -  0(1)  in  L2(R0),  m  -  0,  1,  2 


i 

►  . 

- 

i 

H 

trr 

i 

»  . 

R 

■  ■ 

•€ 

r-- 

< 

►  ,  * 


To  see  that  (6.31)  implies  hm  *  0  note  that,  for  y  >  0,  (5.5)  implies 


(6.32) 


n'  1/2 

yt[y(t2  -  y2)  ]“3/4  h  -(t¥  ■  ■ 7^172- 

u  1 

*  |h(p) |2  dp  • 


dy 


The  equality  of  the  two  integrals  follows  from  the  change  of  variable 


y1/2  y 

P  *  n 


It  is  clear  from  (6.32)  that  (6.31)  implies  hm  *  0  when  m  >  0.  The  case 
m  ■  0  is  equally  simple;  see  (5.8). 

Equations  (6.29)  can  be  written 


h  (P) - 


SOT(/pZ4VV  *mn(p)  * 


These  relations  were  derived  for  functions  h  e  Jf  that  satisfy  (6.16), 
(6.17).  But  sums  of  such  functions  are  dense  in  3C,  by  the  unitarity  of 
<t  and  elementary  properties  of  Xg  .  Hence  the  general  case  follows  from 
(6.29)  by  summing  over  n  and  dropping  the  restriction  (6.17).  Thus 


h~(p) 

m 


(6.33) 


n=0  /p2+y  -y 
r  Mm  Mn 


S+  (/pz+y  -y  )  h"t"(/pz+u  -y  )  y  (p) 
nm  r  *m  Mn  n  v  Mm  Mn'  *mn  v 


l 


P 


Ic  is  clear  from  the  last  expression  that  the  sum  is  finite  for  each 
fixed  p  >  0.  Comparison  of  (6.33)  and  (6.9)  reveals  the  representation 


(6.34)  (§h)m(p)  -  l 

VWm(p)  /p+ynf  yn 

which  is  valid  for  all  h  -  {h  )  €  X„  . 

m  0 

The  known  unitarity  of  S  in  Xg  can  be  used  to  construct  an 
alternative  and  simpler  representation  of  S.  To  derive  it  note  that 
fi+  *  S"1  h~  and  hence  an  analogue  of  (6.34)  for  S~l  *  §*  is  obtained  by 

mm  A. 

interchanging  ^  and  . ,  and  hence  h^  and  h  ^ ,  in  (6.33).  The  result  is 


(6.35) 


(S  g)m(p) 


u  <ar(p) 
Hn 


^p  Z+ym~yn 


8n(v^Vyn) 


It  follows  that  for  all  g,h  e  X0  one  has 


(g,Sh)  -  (S *g,h)  -  [ 
m=*0 


r 


o  (S  g)m(p)  hm(p)  dp 


(6.36) 


Jo  Jo  Xan(P)  SFZTTi 


Snm(’/pi+ym"yn)  gn(/p2+ym"yn)  hm(p)  dp 


m  r  n 


00  fCO 

l 

n=0 


*n<p’> 


l  Xnm(p,)  S^n(p,)  hm(/p  ’ 

m=0 


dp* 


where  the  last  step  follows  from  the  change  of  variable  p2  +  u  *  p’2  +  u 

m  n 


It  follows  that 


(6.37) 


<Sh)n(p)  “  _L  *nm(p)  Sm(p)  hm(^V^) 
m*0 

“  I  s~  (p)  h  (/p^+y  -U  )  . 

“  2  /  \  nm  m  n  m 

VVP> 

The  identities  (6.34)  and  (6.37)  clearly  imply  the  relations 


(6.38) 


The  unitarity  of  S  imposes  further  restrictions  on  the 

coefficients  (S“  }.  The  analogue  of  (6.37)  for  S-1  *  S*  is  obtained 
nm 

from  (6.37)  by  replacing  S  by  S+  .  If  in  (6.37)  h  is  replaced  by 

mn  mn 

S~lh,  represented  in  this  way,  one  finds  after  simplification  using 


7.  Acoustic  Signal  Structure  in  Waveguides. 

The  signals  generated  by  sources  that  are  localized  in  a 
cylindrical  portion  {X  »  (x,y)  :  y0  <  y  <  y0  +  60)  of  the  waveguide  ft 
are  analyzed  in  this  section.  To  make  explicit  the  dependence  on  y  the 
initial  values  will  be  assumed  to  have  the  form 


(7.1) 


u(0,X)  -  f(X,y0)  -  f0(x,y-y0) 


*  g(X,y0)  -  g„ (x,y-y0) 


where 


(7.2) 


supp  f0  U  supp  g0  C  G  X  [O,50 ]  . 


For  simplicity,  the  special  case  for  which  g  *  0  will  be 
discussed.  The  initial  state  (7.1)  with  g  -  0  generates  a  signal  in  the 
simple  waveguide  Q0  whose  asymptotic  momentum  distribution  F0  is  given 
by  (5.14)  with  h  -  f(*,y0);  i.e.. 


(7.3) 


fo(p)  ■  {F0j(p)}  -  {fj(p,y0)> 


where 


(7.4) 


f  j  (p»y0)  -  |  (pJtt.P)  f(X,y0)  dx 


By  (3.1), 


(7-5)  ^°(X,p)  -  (27T^i/2  eipy  <j>j(x)  +  pr-  e~ipy  ^(x) 
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Substituting  this  in  (7.3),  (7.4)  gives 


(7.6) 


Fflj(p>  -  Ffr(P)  +  F^efl(p) 


where 


(7.7) 


Ffr(P) 


Frfl(-p) 


(2tt) 


5^1 


ftn 


e“lpy  ^ (x)  f (X,y0)  dX 


F<?ir(p)  and  FjS^(p)  represent  the  momentum  distributions  generated  by 


f(*,y0)  in  a  doubly  infinite  cylinder  G  x  R  which  have  positive  and 


negative  momenta,  respectively.  Thus,  for  ft0  ,  F?  represents  the 
signal  component  that  propagates  directly  outward  while  Fj  represents 


the  component  that  is  reflected  at  the  end  y  =  0  and  then  propagates 
outward. 

Now  consider  the  signal  generated  in  the  compound  waveguide 
ft  ■  ft0  U  K  by  the  same  initial  state  (7.1)  with  g  =  0.  It  has  an 
asymptotic  momentum  distribution  F  that  is  given  by  (5.36)  with 
h  *  f(*,y0);  i-e.. 


(7.8) 


F (p)  *  (Fj(p)}  *  {f~(p,y0) } 


where 


(7.9) 


fj(p.y0>  *  j 


ft 


^(X,p)  f(X,y0)  dX  . 


Substituting  the  decomposition  (5.19)  into  (7.8),  (7.9)  gives 


Fj (p)  -  Fjir(P>  +  FjC(p) 


(7.10) 
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where  F°  (p)  is  given  by  (7.7)  and  is  again  interpreted  as  the  direct 
part  of  Fj (p) .  F®c(p)  may  be  interpreted  as  the  component  of  F^(p)  that 
is  due  to  scattering  by  the  resonator  K  ■  ft  -  ft0 .  A  short  calculation, 
using  (5.19)  and  (6.37),  gives 


(7.11) 


FjC(P)  -  (S  Frefl)j(p)  +  Fj (p) 


where  F^ (p)  contains  the  exponentially  damped  terms  of  (5.19): 


_  ,  rya-H5  f  *  —  -y/u_-wi(p) 

Fj(p)  -  (2if^  L  *  ♦„<*> 

Jya  JG  (m*0 

(7.12) 


(7.11)  and  (7.12)  suggest  that 


(7.13) 


Fsc  ,  £  Frefl  +  0(1)  in  JC0,  y0  -  »; 


that  is. 


(7.14) 


lim  »FSC  -  S  Frefll-  -  0  . 

y0~°  0 


It  is  not  difficult  to  show  from  (7.12)  that  F^(p)  ■+•  0  when  yg  ■*  ®, 
uniformly  for  p  on  compact  subsets  of  the  set  R  -  '•  m, j *  0,1,2, ••• } 

cf.  [12,  Lemma  14.1].  However,  this  is  not  strong  enough  to  prove 

(7.14) .  Unfortunately,  (7.14)  cannot  be  proved  without  additional 
information  about  ^j(X,p).  This  problem  was  studied  in  [12]  for  the 
analogous  case  of  diffraction  gratings.  There  it  was  shown  [12, 

Theorem  15.1]  that  if  the  analytic  continuation  of  the  resolvent  (A-  z)  1 
across  the  spectrum  of  A  has  no  singularities  at  the  points  z  ■  um  then 

(7.14)  holds.  The  same  method  can  be  applied  to  the  waveguide  problem 
of  this  paper.  The  method  also  extends  to  the  case  of  general  initial 


values  (7.1) ,  (7.2) . 
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